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ABSTRACT 

Merging neutron stars (NSs) are hot candidates for the still enigmatic sources 
of short gamma-ray bursts (GRBs). If the central engines of the huge energy release 
are accreting relic black holes (BHs) of such mergers, it is important to understand 
how the properties of the BH-torus systems, in particular disc masses and mass and 
rotation rate of the compact remnant, are linked to the characterizing parameters of 
the NS binaries. For this purpose we present relativistic smoothed particle hydrody- 
namics simulations with conformally flat approximation of the Einstein field equations 
and a physical, non-zero temperature equation of state. Thick disc formation is high- 
lighted as a dynamical process caused by angular momentum transfer through tidal 
torques during the merging process of asymmetric systems or in the rapidly spinning 
triaxial post-merger object. Our simulations support the possibility that the first well- 
localized short and hard GRBs 050509b, 050709, 050724, 050813 have originated from 
NS merger events and are powered by neutrino- antincutrino annihilation around a relic 
BH-torus system. Using model parameters based on this assumption, we show that 
the measured GRB energies and durations lead to estimates for the accreted masses 
and BH mass accretion rates which are compatible with theoretical expectations. In 
particular, the low energy output and short duration of GRB 050509b set a very strict 
upper limit of less than 100 ms for the time interval after the merging until the merger 
remnant has collapsed to a BH, leaving an accretion torus with a small mass of only 
~ 0.01 Mq. This favors a (nearly) symmetric NS+NS binary with a typical mass as 
progenitor system. 
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1 INTRODUCTION 

Merger events of NS+NS or NS+BH binaries do not only 
belong to the strongest known sources of gravitational wave 
(GW) radiation, they are a lso widely favored as origin of the 
class of short, hard GRBs JBljnnikove^lJ|W8^JPaczvnskJ 
19861: lEichler et ail Il989t IPaczvnskil Il99li iNaravan et al. 
1992). The recent first good localizations of short bursts, 
GRB 050509b, 050709, 050724, 050813, by the Swift and 
Hete satellites at redshifts between 0.160 and 0.722 (see 
iFox et al.ll200l iGehrels et al.ll2005l IVillasenor et al.ll2005l 
and references therein) were interpreted as a possible con- 
firma t ion of this hypothe sis l|Fox_ej 1 j^.| j2005t iBloom et alJ 
120051: iHiorth et all l2005allbl iLee et al] ho05h . because the 
bursts have observational characteristics which are differ- 
ent from those of long GRBs, but which are in agreement 
with expectations for compact object mergers. 

The central engines of such bursts are still poorly un- 
derstood and observationally undetermined. But it seems 



unlikely that the energies required for typical short GRBs 
are set fr ee during the dyna mical phase of the merging of 
two NSs IRuffert et alJll996T) . The production of GRBs by 
the neutrino emitting, hot post-merger NS is also disfavored, 
because the high mass loss rates in a neutrino-driven wind, 
which is caused by neutrino energy deposition near the NS 
surfac e, rules out the production of high Loren tz factor out- 
flow JWooslev fc Baronlll99llWooslevlll993£]) . Instead, the 
long-time accretion of a BH formed from a transiently stable, 
supramassive or hyper massive merger remn ant (for defini- 
tions of these terms, see lMorrison et al. 2004) is a much more 
prom i sing source (e.g.. IWooslevI Il993bl IRuffert fc Janka 
I 1999L iPopham et"all ll999L iRosswog etal] l2003l ILee et al. 
2005), provided the BH is surrounded by a sufficiently mas- 
sive accretion torus. Due to the geometry of the BH-torus 
system with relatively baryon-poor regions along the rota- 
tion axis, thermal energy release preferentially above the 
poles of the BH by the annihi l ation of neutrino-antineu trino 
iyv) pairs ll Jaroszvriskil Il99"3. iMochkovitch et al]ll993^l can 
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lead to collimated, highly relativistic jets of baryonic mat- 
ter with properties in agreement wit h those needed to ex- 
plain short GRBs (lAlov et alJ I2005T) . Ultrarelativistic jets 
were found to develop if the rate of thermal energy deposi- 
tion per unit solid angle is sufficiently large. Alternatively 
or in addition, magneto-hydrodynamic (MHD) field amplifi- 
cation in the rapidly spinning disc could drive polar jets, or 
magnetic coupling between the rotating BH and the girding 
accretion torus could tap the rotational energy of the BH 
felandford fc Zn aiek||l977| ) and could help powering MHD- 
drive n outflow (e.g. jBrown et al.l200d . lDrenkhahn fc Spruid 
120021) . 

Torus mass as well as BH mass and rotation are thus 
crucial parameters that determine the energy release from 
the merger remnant on the secular timescale of the vis- 
cous evolution of the BH-torus system. In this letter we 
present a first set of results from new three-dimensional 
(3D), relativistic smoothed particle hydrodynamics simu- 
lations of NS mergers with a physical equation of state 
(EoS), which aim at establishing the link between the rem- 
nant properties and t hose of the binary systems . Previ- 
ous such attempts fsee lBaumgarte fc Shapiral2003l for a re- 
view) were either per formed in the Newton ian limit (e.g., 
[Ruffert fc Jankal200fl and references therein; iRosswog et alJ 
2003), or with polytropic or otherwise radically simplified 
treatments of the high- and low-dens ity EoSs and their 
temperature-depend ent behaviour fe.g.. lShibata et al]|2003l 
IShibata et all [20051. or they considered disc formation by 
viscous or magnet ic effects during the secular evolut ion of 
hypermassive NSs (iDuez et all2004lDuez et all21)l )■',:. Here 
we instead highlight torus formation in NS mergers as a dy- 
namical process and argue that this is consistent with the 
recent observations of short GRBs. 



2 NUMERICAL METHODS AND MODELS 

We employ an improved version of our relativistic smoothed 
particle hydrodynam ics (SPH) code llOechslin et al.l I2002L 
lOechslin et ai]l2004) . which solves the relativistic hydrody- 
namics equations together with the Ein stein field equation 
in the conformally flat approximation dlsenberg fc Nesterl 
Il98d. IWilson et al.lll996h . The code now allows for the use 
of a tabulated, non-zero temperature EoS and solves the en- 
ergy equation in a form without expli cit time derivatives of 
the metric elements on the RHS (cf. lOechslin et al.l l2002t 
eqn. A8). The properties of the stellar plasma are describ ed 
by the non-zero temperature EoS of IShen et alJ Il998al lbl 
which is used with baryonic density, internal energy and 
electron fraction (proton-to-baryon ratio) Y e as input, and 
pressure and temperature as output. Since the backreaction 
of the neutrino emission on the stellar fluid is small on the 
timescales considered here, we ignore neutrinos in our mod- 
els. As a consequence, the electron fraction Y e needed by the 
EoS as an input remains constant on the Lagrangian SPH 
p articles. Details on the numerical scheme will be presented 
in lOechslin et alJ J2006ft 

We start our simulations shortly before the tidal insta- 
bility sets in and follow the evolution with typically 400'000 
SPH particles through merging and torus formation until 
either the collapse of the merger remnant sets in, or a quasi- 
stationary state has formed. 



Initial conditions for our simulations are generated by 
placing two NSs in hydrostatic equilibrium on a circular or- 
bit of a given orbital distance and by relaxing the configu- 
ration with the use of a small damping force into a circular 
orbital motion with an irrotational spin state. The orbital 
velocity is adjusted during this process in order to obtain 
the binary in orbital equilibrium. The initial electron frac- 
tion F e is obtained by requiring neutrinoless /3-equilibrium 
for cold neutron star matter, i.e., we determine Y e from the 
condition that the (electron) neutrino chemical potential 
is equal to zero for a given density (and initial temperature 
T ~ 0). 

We consider a variety of models with different NS 
masses and mass ratios. Other unknown parameters like the 
NS spins and the EoS are kept fixed except for two mod- 
els where thermal effects in the EoS are neglected. This is 
realized by reducing the EoS to the twodimensional slice at 
T — with density and Y e as input and pressure and in- 
ternal energy as output. This reduction has no influence as 
long as shocks are absent, i.e. as long as the fluid evolves 
adiabatically. In the presence of shocks, however, the T = 
case corresponds to the extreme situation that a very effi- 
cient cooling mechanism extracts immediately the entropy 
and internal energy generated in shocks. 

In Table the key parameters characterizing our differ- 
ent models are summarized. We have investigated NS+NS 
binaries with mass ratios q between 0.55 and unity, roughly 
covering the theoretical ra nge obtained by st ellar popula- 
tion synthesis calculations feulik et al]l2003h . Predictions 
of the mass distribution, however, are sensitive to parame- 
ters that govern single and binary star evolution and still 
contain significant uncertainties. A comparision with obser- 
vations of galactic double NS binaries, which currently span 
a range of q- val ues between 0.67 and nearly uni ty (see the re - 
cent reviews by lLattimer fc PrakashlEool and IStairdl200^ , 
is therefore difficult and is also hampered by the still small 
number of detected NS+NS systems. 

In the following analysis and discussion of our models, 
we have synchronized the time axis to the time of maxi- 
mal gravitational wave luminosity. This allows for a better 
comparison of the temporal evolution. 



3 RESULTS OF MERGER MODELS 
3.1 Dynamics of Merging 

The evolution of a NS binary during the adiabatic inspiral is 
driven by gravitational radiation reaction. Once the binary 
becomes unstable to tidal forces, the hydrodynamic evolu- 
tion sets in, leading to the tidal stretching of one or of both 
companions and to the formation of a central hypermas- 
sive NS surrounded by a thick, neutron-rich accretion torus. 
Depending on its mass and its rotation rate and the EoS, 
the hypermassive NS will collapse immedi ately or after a 
(shor t?) delay to a rapidly rotating BH (e.g. lMorrison et alJ 
12004 . 

The dynamics and post-merging structure largely de- 
pend on the mass ratio q of the binary. In asymmetric sys- 
tems with q significantly smaller than 1, the less massive but 
slightly larger star is tidally disrupted and deformed into an 
elongated primary spiral arm, which is mostly accreted onto 
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Table 1. Parameters of the considered NS+NS models. M\ and M2 denote the individual gravitational masses of the NSs in isolation, 
while Msum = Mi + M2 stands for the sum of the two. Note that the total gravitational mass M is slightly smaller than M SU m because 
M also involves the negative gravitational binding energy between the two stars. Mo is the total baryonic mass, q = M1/M2 and 
qm = Mp. 1 /Ma. 2 a re the gravitational and baryonic mass ratios, respectively. 'Shen' stands for the full, non-zero temperature EoS of 
IShen et al.ll998albl while 'Shen_c' denotes the T = slice of this EoS table. Characteristic data of the post- merger system are read off at 
t ~ 6ms. Mo i( jisc arl d -/disc ar e the baryonic mass and the angular momentum of the disc, respectively. M Tern is the remnant gravitational 
mass, and Orena = (-/total — Jdisc)/Mf Bm is the corresponding spin parameter, jisco an d ISCO are the specific angular momentum and 
the radius of the ISCO of a Kerr-BH with the gravitational mass and spin parameter of the remnant. 
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the more massive companion. Its tail, however, contributes a 
major fraction to the subsequently forming thick disc/torus 
around a highly deformed and oscillating central remnant 
(see Fig. HjJ panels (a) and (b)). In systems with q ~ 1 and 
nearly equally sized stars, both stars are tidally stretched 
but not disrupted and directly plunge together into a de- 
formed merger remnant (see Fig. . In all models, excited 
radial and non-radial oscillations of the compact remnant 
periodically lead to mass shedding off the surface and to 
ejection of material into secondary tidal tails during the 
subsequent evolution. As illustrated in Figs. |S| and [7| this 
effect takes place both in symmetric and asymmetric cases, 
but is stronger in the latter because of the larger triaxial 
deformation of the post-merger object. 

3.2 Disc Formation 

If and when the compact post-merger object collapses to a 
BH (the secular evolution driven by GW radiation reaction, 
viscosity and MHD cannot be followed by 3D simulations), 
pressure support from the central remnant to the disc will 
vanish and matter will be prevented from infall mostly by 
rotational support. As a consequence, we distinguish future 
disc matter from the central remnant by demanding the 
specific angular momentum j to be larger than the one 
associated with the innermost stable circular orbit (ISCO) 
of a Kerr BH with the gravitational mass and the spin 
parameter of the central remnant. 

The ISCO can be analytically determined in the ca se of 
the Boyer-Lindquist Kerr metric feardeen et al]ll972lh In 
our case we shall use an approxima tive pseudo-Kerr metric 
(see, e.g.. lOrandclement et alJl2002t) which is both isotropic 
and conformally flat, consistent with the coordinates we used 
for our numerical simulations. The ISCO can then be found 
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Figure 1. Disc formation in models with non-zero temperature 
EoS and Af sum = 2.8 Mq. Plotted is the disc mass as inferred 
from our disc mass criterion (see text). The final disc mass as 
given in Table HI is read off at about 6ms. Clearly visible is the 
rapid rise at t ~ 0, which depends strongly on the mass ratio. On 
the other hand, the further evolution after merging is similar in 
all cases. 



among all circular orbits in the orbital plane by minimizing 
the specific angular momentum along the radial coordinate. 

Using this criterion, we are now able to identify the 
disc matter and disc mass for the various models consid- 
ered. Their very different merger dynamics is reflected in 
the different evolution of the disc mass (see Fig. [TJ. The 
rapid rise in the asymmetric models around t ~ 0ms can be 
associated with the development of the primary spiral arm. 
After this dynamical phase shortly after merging, a nearly 
stationary torus forms and the disc mass settles down to a 
slowly changing value, which is plotted versus mass ratio q 
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Figure 2. Disc masses versus mass ratio q for all non-zero tem- 
perature runs. 
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Figure 3. Same as in 

Fig. EI but comparing the differences be- 
tween T = EoS and non-zero temperature EoS for two different 
mass ratios q. 



for all models in Fig. [5] We find a rapid rise of the disc mass 
when q drops below unity, and a flattening for q-values be- 
low ~ 0.8. The disc mass increases from about O.O5M0 at 
q = 1 to about O.26M at q = 0.55. 

We checked the consistency of our disc mass determina- 
tion by comparing the obtained masses with the amount of 
matter residing outside the equatorial radius of the remnant 
r rem . We determine r rcm in a radial density profile by locat- 
ing the point between the steep density gradient at the rem- 
nant surface and the slower fall-off in the disc. In all cases, 
we found A/ gas (r > r rcm ) > M gas (j > jrsco) = Mo.disc, 
suggesting that gas pressure still plays an important role 
in the quasi-stationary phase. Our disc masses are there- 
fore lower limits, since pressure support will not completely 
vanish even after BH formation. 

In Fig. HJwe have color-coded for a T = and 
model in red the material that belongs to the future disc, 
whereas the material that currently fulfills the disc crite- 
rion is shown in yellow. By definition, these two attributes 
coincide at the end of the simulation. We see that the fu- 
ture disc matter originates mainly from the spiral arm tips 



simulation 

postprocessing 




Figure 4. Total angular momentum of a blob of matter in the 
primary spiral arm of Model C1216. The dashed line represents 
the results of the simulation, the solid line those of a postprocess- 
ing analysis. 



and from the stellar surface. This indicates that the forma- 
tion of a post-merger disc is linked to the presence of spiral 
arms. Indeed, the yellow tips of these tidal tails in Fig.|S|and 
also the enhanced specific angular momentum in the outer 
edges of the tidal tails in Fig. ^suggest angular momentum 
transfer from the central remnant to the spiral arm tips via 
tidal torques. This hypothesis is confirmed by a postprocess- 
ing analysis of Model C1216 where we integrated the torque 
exerted by gravitational and pressure forces from the com- 
pact deformed remnant on the yellow mass elements in the 
primary spiral arm as a function of time. The result for the 
total angular momentum evolution agrees very well with the 
simulation result (Fig.^J. 

An interesting effect can also be seen by comparing 
models with T = and T / EoSs (Figs. 01 U E3- In 
Fig. [3] we show the disc mass evolution of the cold Models 
C1216 and C1315 together with their non-zero temperature 
counterparts S1216 and S1315. Obviously, the first increase 
in disc mass, which can be associated with the cold primary 
spiral arm, is unaffected by temperature effects. Thermal 
pressure effects during the following evolution, however, re- 
duce the further increase in disc mass in the non-zero tem- 
perature case, whereas in the cold case the disc mass in- 
creases by an additional ~ O.O7M0. We determine two pos- 
sible reasons for this difference. On the one hand, the core of 
the cold merger remnant contracts more strongly due to the 
absence of thermal pressure. Therefore, more gravitational 
binding energy is converted into rotational energy and the 
rotation rates become higher, supporting the formation of 
spiral arms. On the other hand, the shock-heated matter 
from the collision region between the two NSs forms a halo 
of hot matter which engulfs the compact remnant-torus sys- 
tem in the T / case (Figs. 03 EJ. This halo may damp 
the oscillations of the remnant and the development of large 
deviations from axisymmetry. 
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Figure 5. Density p (left), specific an- 
gular momentum j (middle) and temper- 
ature distribution T (right) in the orbital 
plane for the T / Model S1216 (top) 
and the T = Model C1216 (bottom, 
without temperature). The black dotted 
line in all panels and the black solid line 
in the color bars for j correspond to the 
location where the specific angular mo- 
mentum is equal to the value of Jisco at 
the given time (see text). Note that the 
color bar for the upper right panel is lim- 
ited to 10 MeV such that the temperature 
distribution in the disk becomes visible. 
The maximum temperature in the merger 
remnant rises up to ~30 MeV at the given 
time. 



4 CONSTRAINTS FROM OBSERVED GRBS 

In this section we shall discuss how the properties of 
the GRB engine may be constrained from the recent ob- 
servations of four short-duration GRBs. To this end we 
shall make the assumption that the bursts originated 
from ultrarelativistic jets whose acceleration was mainly 
driven by thermal energy, provided by the annihilation of 
neutrino-antineutrino iyv) pairs in the vicinity of post- 
merger BH-torus sy stems (see, e.g., iRuffert fc Jankalll999l 
Scti awan et al . 2004). Although we will refer to this scenario 
in our discussion because numerical models provide at least 
order of magnitude information of involved parameters, our 
arguments are more general and apply basically also to mod- 
els that consider MHD-powered GRB jets. 



4.1 BH- Torus System Parameters 

Linking the measured apparent (i.e., isotropic-equivalent) 
energy -E 7l i so of a short GRB to the energy output from the 
central engine and the mass A/ acc accreted by the BH (dur- 
ing the phase when the neutrino emission from the accre- 
tion torus is sufficiently powerful to drive the jets) involves 
a chain of efficiency parameters corresponding to the differ- 
ent steps of physical processes between the energy release 
near the BH and the gamma-ray emission at ~ 10 14 cm : 



-E 7 ,iso = fi h h fn 1 U M^c 



(1) 



Here /i denotes the efficiency at which accreted rest mass en- 
ergy can be converted to neutrino emission, fi is the conver- 
sion efficiency of neutrinos and antineutrinos by annihilation 
to e pairs, fa is the fraction of the e ± -photon fireball en- 
ergy which drives the ultrarelativistic outflow with Lorentz 
factors T > 100 as required by GRBs, fn = 2fij e t/(47r) = 
1 — cos 0j e t denotes the jet collimation factor defined as the 
fraction of the sky covered by the two polar jets (with semi- 
opening angles f?j ct and solid angles fij e t), and fi is the frac- 
tion of the energy of ultrarelativistic jet matter which can be 
emitted in gamma rays in course of dissipative processes that 
occur in shocks when optically thin conditions are reached. 

These parameters are constrained to some degree by nu- 
merical and analytic work, but their values are still rather 
uncertain and might vary strongly with time-dependent 
and system-dependent conditions of the source. "Typi- 
cal" values from merger and accretion simulatio ns are: 
fi - 0.05 jLee et alj 120051 [Setiawan et alJ 120041) f 2 ~ 
0.001 ... 0.01 iRuffert fc Jankalll999l. ISetiawan et all 120041) . 
/a ~ 0.1, fn - 0.01 ...0.05 and .f 4 <0.2 from est i mates 
for internal shock models foa igne fc Mochkovitchl 1 19981 . 
iKobavashi fc Sad 1200 ll . iGuetta et alJl200ll and references 
therein). 

The cited values for /a and fn were determined by rela- 
tivistic hydrodynamical calculations of the formation and 
propagation of jets driven by thermal ene rgy deposition 
around BH-torus systems llAlov et al.ll2005l) . These simu- 
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t=9-673ma t=9.6655ms 

(e) (f) 

Figure 6. Characteristic evolutionary phases of Model C1216 (left panels) and Model S1216 (right panels). Panels (a) and (b) show 
the merging phase with a primary spiral arm forming. Plotted is every 10th SPH particle around the equatorial plane. The matter of 
the two stars is represented by green and blue particles, respectively. Particles that end up in the disc at the end of the simulation are 
marked in red and the ones that currently fulfill the disc criterion (see Sect. 13.21 are plotted in yellow. Panels (c) and (d) display the 
situation ~ 2 ms after merging, when a secondary spiral arm appears. Panels (e) and (f) show the time when a quasi-stationary torus 
has formed. Note that the particle density is not strictly proportional to the density of the fluid because the particles near the initial 
stellar boundaries are attributed a smaller particle mass. These particles preferentially end up in the torus so that the particle density 
there appears enlarged. 
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lations revealed that about 10-30% of the e ± -pair plasma 
fireball energy are used for hydrodynamically accelerating 
the baryonic matter in the jets to ultrarelativistic velocities. 
Typical jet half-opening angles were found to be between 
10° and 15° corresponding to fn ~ 1.5 x 10" 2 ...3.4 x 10 -2 
l|Alov et al]l200SL IJanka et al J (20051) . Two of the four well- 
localized short G RBs indeed prov i de hints for this d egree 
of collimation (see lFox et aljEooil iBerger et alJl2005h pre- 
dicted by the theoretical work. 

Because of the lack of detailed information about how 
the factors f\ to fi and fn depend on the properties of 
NS+NS binaries and their relic BH-torus systems, we made 
the simplifying and bold assumption that the set of param- 
eter values is the same for all cases, (/i, fe, fz, fn, fi) = 
(0.1, 0.01, 0.1, 0.01, 0.1), and applied Eq. JTJ to the four ob- 
served short GRBs, using the measured isotropic-equivalent 
gamma-energies of Table H Doing so, we obtained the es- 
timates listed in that table for the accreted mass Macc 1 . 
Assuming further that the duration i 7 of the GRB (at the 
location of the redshifted source) is an upper bound for 
the duration of the a ccretion of A/ acc by the BH, t acc < i 7 
(cf. lAlov et alJliooBT) . we can deduce lower limits for the 
average mass accretion rates, 

Mace > ^ • (2) 

These are also listed in Table |2] and plotted versus M acc 
m Fig. G3 We obt ain in all cases values for the torus 
masses in the range of those determined in this work 
(Fig. and values for the mass accretion rates in the 
ballpark of theoretical expectations for post-merger BH 
accretion. Simulations of the latter find mass accretion 
rates of typic ally between several O . ltfas" 1 and several 
Mp, s" 1 (e.g.. IRuffert fe Jankalll999l ISetiawan et alJ 12004 
iLee fe Ramirez-Ruijl200l iLee et alll2005l) . Of course, the 
numerical values for M acc and M aC c in Tableware very un- 
certain and have to be taken with caution. Shifts by fac- 
tors of a few are easily possible because of probable system- 
dependent variations of the efficiency factors /i to fi and 
of the jet collimation factor fn- Nevertheless, it is very en- 
couraging that the simple and straightforward calculations 
based on Eqs. Q and lead to reasonable numbers in all 
cases. 

GRB 050509b, for example, is particularly interesting 
because of its very short intrinsic duration of t~, ~ 33 ms 
and an extraordinarily low isotropic energy of a few 10 4S erg 
if it is located at z = 0.225 iGehrels et al.l2005h . Such a low- 
energy event can be explained by the accretion of discs in the 
lower range of masses found here for symmetric NS+NS bi- 
naries, even when rather conservative assumptions are made 



1 Note that the (initial) torus mass may be larger than the ac- 
creted mass that we determine from the observed GRB proper- 
ties. A part of the torus may fall into the BH without producing 
a release of energy which is strong enough to power ultrarelativis- 
tic GRB-jets. Some part of the torus mass will also be ejected by 
wind s driven by neutrino energy depo sition (see e.g.. IRuffert et aO 
Il997l IRosswog fe Liebendorfeij|200^. v iscous heating, or MHD 
effects (see, e.g., iDaigne fe Mochkovitchll2002D . and some part of 
the torus mass may be lost equatorially due to the outward trans- 
port of angular momentum by viscous shear in the (magnetized) 
accretion torus. 



for the incompletely known efficiencies of the various pro- 
cess steps between the energy release near the BH and the 
gamma-ray emission. We find that an accreted mass of 10 -3 - 
10~ 2 Mq is well sufficient to account for the energetics of 
GRB 050509b. Of course, the mass could also be larger if 
the energy conversion to ultrarelativistic outflow and GRB 
emission is less efficient than we assumed. 



4.2 Constraints from the low-energy, 
short-duration GRB 050509b 

If further observations substantiate the link between NS+NS 
mergers and short GRBs, this may set constraints on the 
post-merging evolution and the nuclear EoS. On the one 
hand, the hypermassive NS should escape the collapse to 
a BH for a sufficiently long time so that an accretion torus 
can form by the described dynamical processes. On the other 
hand, the collapse of the hot, neutrino radiating NS should 
not be delayed too much, because neutrino energy depo- 
sition in the surface-near layers of the N S will lead to a 
massive baryonic wind (Duncan et al. 1986. | Wooslevlll993al . 
iQian fe Wooslevll996l.fThompson et al 120011b which can se- 
riously endanger the subsequent formation of a GRB jet or 
fireball. 

The mass loss rate of the neutrino radiating neutron 
star due to this nonrelativistically expanding outfl ow of 
bayonic matter was e stimated by iWooslevI [l993a) and 
jQian fe Wooslevlll996ri and is approximately given by (see 
iHeger et alJl200oT) 

M wind » 8.5 x 10" 3 L^p, i6 a R% 3 M 3 " 2 M s" 1 , (3) 

where !/„,,+[?,,, 53 is the luminosity of electron neutrinos plus 
antineutrinos in units of 10 53 ergs -1 , R30 the NS radius nor- 
malized to 30 km, and M3 the NS mass in units of 3 Mq . 
This mass loss rate can be up to ~ 10 -2 Mq s _1 for the 
neutrino luminosities found in NS+NS merger simulations 
and the masse s and radii of the compact post-merger ob- 
ject ( see IRuffert fc Jan kal l200lL IRosswog fc Ramirez-Ruizl 
20021, IRosswog fc Liebendorferll2003D . Upon colliding with 
each other, the two merging neutron stars heat up by 
shocks and compression to central temperatures of sev- 
eral 10 MeV (Figs. |5] |5J . Within only a few millisec- 
onds after the final plunge, the neutrino emission reaches 
luminosities of 10 53 ergss~ 1 or h igher (for more infor- 
matio n, see IRuffert et all Il997t IRuffert fc Jankal l200l| 
I Rosswog fc Liebend5rfe? 120031 . iRosswo^^^R^mirez^^izl 
2003), and the post-merger object should start losing mass 
in the neutrino-driven wind 2 . Equation however, does 
not take into account the effects of the very rapid rotation 
of the neutron star on the wind. This question is still un- 
explored, but the rapid spinning of the remnant born in a 
NS+NS merger event may lead to changes of the functional 

2 The development of the neutrino-driven wind is not seen in 
current NS+NS merger simulations, because it requires a suit- 
able treatment of neutrino transport, whereas the present models 
describe neutrino effects with a trapping scheme at best, which 
releases neutrinos from the stellar medium according to their lo- 
cal production rates, properly red uced due to a finite diffusion 
times cale (see IRuffert et alj Il996l and IRosswog fe Liebendorfeil 
120031) . 
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(a) 



(b) 



Figure 7. Same as panels (a) & (c) and (b) & (d), respectively, of Fig. El but for Model S1414. The primary spiral arm is absent in the 
symmetric binary case and the post-merging (secondary) arms are smaller than in the asymmetric models 



(a) 



T=0 15 30MeV 










T = -1.40064 ms 


T=0 15 30MeV 










T = 0.600967 ms 




(c) 

Figure 8. Four stages of the NS+NS merger Model S1216. The surface is chosen to correspond to a density of 10 10 gcm -3 , the 
temperature distribution is visible color coded in the octant cut out from the three-dimensional mass distribution. Temperatures up to 
30 MeV are reached at the center only milliseconds after the two neutron stars have merged. Time is given in the lower left corner of the 
panels, and length is measured in units of 1.47km (top right corner of each panel). 
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Figure 9. Estimated accreted masses, M, 
for the average mass accretion rates, M, 
localized short GRBs of Tabic |2] 



Table 2. Estimated accreted masses, Af acc , and lower bounds 
for the average mass accretion rates, M acc , assuming post- merger 
BH-torus systems being the energy sources of the recently ob- 
served, well-localized short GRBs. The estimates are based on 
Eqs. Q and J2J- The quantity i? 7i i so is the isotropic-equivalent 
7-ray burst energy (corrected for the cosmological redshift z of 
the burst), and t~ is the GRB duration at the source, computed 
as t-y = Tgo/(l + z) from the measured 90%-inclusive interval 
of hig h-energy e mission. The observational data were taken from 
iFox et alJ <2005lV 
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dependence of M on the NS properties and may also imply 
a strong pole-to-equator difference of the mass-loss rate. 

The nonrelativistic wind will create an extended bary- 
onic halo that surrounds the merger remnant, before the 
compact central part undergoes collapse to a BH and the 
outer, high-angular momentum matter assembles in a torus 
around the BH. A jet launched from such a system and ex- 
panding into the pre-existing wind halo may then sweep up 
so much mass - even if the halo has a low density and small 
mass - that the high Lorentz factor outflow needed for GRBs 
(r>100) is prev ented. This was, f or example, seen in the 
type- A models of lAlov et, alJ (|2005ll . A rather fast collapse 
of the relic hot, supramassive or hypermassive remnant to 
a BH is therefore required, if the cloud of baryonic matter 
produced by the neutrino-driven wind from the transiently 
stable NS shall not become such a hazard for the GRB 3 . 

3 Note that neutrino energy deposition will lateron also drive 
a baryonic wind off the accretion torus, after the latter 
was heated by viscous dissipation of rotational energy and 



A corresponding mass limit for the wind halo around 
the newly formed BH-torus system can be derived by the 
following arguments. In order for the baryon loading to re- 
main sufficiently low so that acceleration to Lorentz factors 
above a value To (100 or more) is possible, the mass of the 
halo swept up by two relativistic GRB jets, M awep t, is con- 
strained by the condition I?km(r > ro)/(M awe ptC 2 ) > IV 
Here i?kin(r > To) is the kinetic energy of matter which 
can accelerate to Lorentz factors larger than To in the two 
polar jets that expand away from the BH in the axial di- 
rection. For a jet collimation factor /n, this kinetic energy 
is linked to the isotropic equivalent 7-energy of the GRB 
by Skin = -E T ,iso/n/ fi, when fi is again the efficiency for 
the conversion of kinetic energy to gamma rays. On the 
other hand, the wind mass swept-up by the jets is given 



as M E 



wept 



= /ftaM, 



where M w 



Af„i n dtNS is the 



mass of the wind for a lifetime £ns of the neutron star, and 
a < 1 is a fudge factor which accounts for the fact that 
the wind from a rapidly rotating neutron star might be less 
dense along the rotation axis than calculated from Eq. ©. 
Combining all, we get the condition 



T f 4 a A/ W i n d c 2 



(4) 



With To = 100, / 4 ~ 0.1, Afwind = 10" 3 M Q s _1 (this as- 
sumption is on the low side of reasonable numbers from 
Eq.OJ, and £ 7 , iso from Tabled we find for GRB 050509b 
(the lowest-energy case and thus the one providing the 
strongest limits) that 4ns < 1ms, assuming a = 1. This 
means that due to the low energy of GRB 050509b even a 
tiny amount of baryonic pollution by swept-up wind mat- 
ter would have prevented ultrarelativistic motion of the jets. 
The neutrino radiating merger remnant therefore must have 
collapsed to a BH essentially immediately after the binary 
NS merging and the onset of the neutrino-driven wind. Since 
it cannot be excluded that the density of the neutrino-driven 
wind along the rotation axis is much lower than estimated by 
Eq. 10, a more conservative estimate with a ~ 0.01 yields 
t N s < 100 ms. 

If GRB 050509b originated from a NS+NS merger, this 
result therefore suggests that the merger remnant was a 
hypermassive object, which was transiently stabilized b y 
its rapid and differential rotation feaumgarte et all hoOO). 
but had a very short lifetime. Referring to the analysis by 
iMorrison et alJ J2004T) - see their Table 2 - this can either 
mean (i) that, if the high-density EoS was stiff, the system 



has started to ra diate neutrinos with high luminosities (see 
the discussion in |Ruffgrt_et^L [119971 I Ruffert &_Jan k3 Il999f 



iRosswog &: Rarnirez^^rrzl ' I^OO? ' R^sswo^^t^lTl2003T ^nd find 
inform ation about the neu trino emission from such accretion 
tori in ISetiawan et alj|2004j and iLee et aljl2005l and references 
therein). This wind does not provide such a hazard for GRB- 
viablc outflows, because it originates from outside the innermost 
stable circular orbit, and angular momentum conservation pre- 
vents it fro m filling the axial r egion s above the poles of the 
BH. In fact, Ros swog et alJ l2003h and lRosswog &: Ramirez- Ruizl 
(2003) consider the pressure of such a "wind envelope" as 
helpfu l for the collimation of a i/p-annihilation-driven axial jet 
(see iLevinson fc Eichlerl EooTj) . an effect which turned out to 
be secondary for getting collimated, highly-relativistic outflow 
from BH-torus sy stems in recent hydrodynamical jet simulations 
iAlov et alJl2005lV 
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mass of the merging binary was unusually high compared to 
the known galactic double neutron stars, which all possess 
combine d masses in the range of M sum = Mi + M2 = 2.6- 
2.8 Mq iStairsl 12004) : or, alternatively, (ii) that the nuclear 
equation of state is so soft that it could not support the 
merger remnant despite of its rapid differential rotation and 
thermal pressure; or (iii) that a very efficient mechanism 
was at work which destroyed the differential rotation of 
the compact remnant and extracted a sizable amount of its 
angular momentum on the required short timescale. This 
could, for example, happen by hydrodynamical effects and 
mass shedding and /o r by strong gravitational wave emission 
JShibata et alJl2005l) . 

Case (i) might be rejected as unlikely. C ase (ii) tends to 
favor the "softest" of the EoSs surveyed bv iMorrison et all 
i2004l) which support nonrotating NSs only up to a grav- 
itational mass around 1.65 Mq. Such EoSs, however, are 
excluded by the recent accurate measurement of a mass 
of 2.1 Mp) for the millisecond pulsar PSR J0751+1807 
jNice et al.ll2005l) and by other neutr on stars with masses 
beyon d 2 Mq (see, e.g., the review bv lLattimer fe Prakasbl 
2004). An EoS which is sufficiently stiff to fulfill this observa- 
tional constraint can also ensure the transient stability of the 
merger remnant of binary NS systems with a canonical mass 
of 2.6-2.8 Mq (see the discussion of the critical mass for hy- 
permassive NSs bv lShibata et al] l2005). Therefore case (iii) 
seems to offer a plausible scenario for GRB 050509b. The 
NS+NS binary had a typical mass and the merger remnant 
was well above the maximum mass of stable, rigidly rotat- 
ing neutron stars. The redistribution of angular momentum 
by non-axisymmet ric hydrodynamic in teraction and gravi- 
tational radiat ion iShibata et all l^OOS 5 ) or by viscosity and 
MHD effects ijShapfrol 1200'fF has then driven the hyper- 
massive object to gravitational instability on a timescale of 
<C 100 ms. The short accretion timescale (i aC c < t-y) suggests 
that the remaining torus mass was small. This is consistent 
with our conclusion that little mass was accreted by the BH, 
which we have drawn from the low energy of GRB 050509b 
(Table |5J . And it is consistent with the requirement coming 
from optical limits th at little radiating material was ejected 
jHiorth et alJl2005al) . e.g. by neutrino-driven winds or by 
mass shedding due to viscous transport of angular momen- 
tum within the accretion torus. In view of the torus forma- 
tion models presented here, this points to a symmetric or 
nearly symmetric NS+NS binary with a typical mass as the 
source of GRB 050509b. 



5 DISCUSSION 

Although our analysis of the four well-localized short GRBs 
was based on a number of unsettled assumptions, we found 
that the measured GRB energies and durations lead to es- 
timates for M acc and M aC c that are roughly consistent with 
expectations from theoretical models of NS merging and 
post-merging evolution. However, in spite of this amazing 
result we think that it would be premature to claim that 
relic BH-torus systems from binary NS mergers and their 
associated neutrino emission can power all short GRBs. 

Short GRBs have durations between a few milliseconds 
and about two seconds, and they might differ in their en- 
ergy output by a factor of 100 or more. Observations also 



reveal large differences with respect to their lightcurves and 
spectral hardness ratios. This wide range of short GRB prop- 
erties leaves plenty of room for the possibility that different 
compact binaries contribute to the observed events. For ex- 
ample, not only NS+NS but also NS+BH mergers can lead 
to BH-torus systems with sizable torus masses and favor- 
able properties for GRB production. Even the physical pro- 
cesses which play a dominant role in the energy release of 
the central engine and jet formation might differ between 
the events. In fact, it might turn out that this is needed 
to account for the observed diversity of short GRBs and 
their afterglows. Flare-like events followed by a very rapid 
decay were, for example, observed in case of GRB 050724 
hundreds of seconds after the initial GRB. This was inter- 
preted as a consequence of long-time activity of the cen- 
tral engine, which might point to a partial disruption an d 
gradual accretion of a NS by a BH fearthelmv et al1l2005l) . 
Episodic and long-lasting mass transfer from the NS to the 
BH was found to be possible for certain NS /BH mass ratios 



in combination with special pro 
state (see e.g.,|Kluzniak & Led 


perties of the NS equation of 
1998l.lPorteeies ZwartJll998l. 


|Lee fc Kluzniakl Il999t Ijanka et alJ Il999l 


iRosswog et alJ 



12004 lRosswog|l2005l . iDavies et al]l2005 

The rather uncertain estimates of the merger rates also 
allow for the possibility that not all NS+NS / BH mergers 
produce GRBs (see, e.g.. lGuetta fc Piran |20 05|. taking into 



Alov et alJ J2OO 



account that the hydrodynamic models by[, 
yield jet semi-opening angles of ~ 10° instead of a value of 
~ 1° as found by Guetta & Piran with the assumption that 

all mergers make GRBs). 

We therefore refrain (unlike othe rs, see, e.g.. iLee et alJ 
120051 IRosswoe fc Ramirez-Rurzll2003l) from making predic- 
tions of the distribution of GRB properties on grounds of 
current simulations of NS+NS or NS+BH mergers. Such 
attempts are seriously hampered by the large number of 
unknowns in the theoretical picture and in particular their 
relation to the properties of compact object mergers: 

• The distributions of binary parameters (NS and BH 
masses and spins, binary mass ratios) are highly uncertain. 
Theoretically, because of our incomplete understanding of 
the binary formation and evolution and of the supernova ex- 
plosions which termina te the lives of ma ssive stars and give 
birth to neutron stars feulik et al.ll20o3) . And observation- 
ally, because of the rather limited sample of known NS+NS 
or N S+BH binaries (e.g.. lLattimer fc Prakash1l2004 Istairsl 
120041) . The latter aspect is particularly bothersome in view 
of the fact that only a minor fraction (of order 10%?) of 
NS+NS/BH mergers might be able to produce GRBs and 
might be needed to account for the observed short bursts. 

• NS+NS/BH merger models require much more work. 
Relativistic simulations of compact object mergers are 



4 Simulations and analytic studies of NS+BH mergers have so 
far been carried out only in Newtonian gravity and disregarding 
the probable rotation of the BH. A Newtonian treatment, how- 
ever, does not reproduce the correct mass-radius relation of a 
relativistic neutron star, for example. Current predictions of the 
merger evolution and of the dynamics of disc formation and ac- 
cretion in such extremely relativistic systems can therefore not be 
considered as reliable and conclusive, neither quantitatively nor 
qualitatively. 
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needed for different (non-zero temperature) nuclear equa- 
tions of state and different parameters of the binary sys- 
tems. A softer nuclear EoS, for example, might lead to a 
relation between binary parameters and torus masses that 
differs from the resu lts of this work, w hich are based on 
the rather stiff EoS of lShen et all lll998allb|) . NS+NS merger 
models must be consistently evolved from the last stages of 
the progenitor system, through the moment of BH forma- 
tion, to the subsequent accretion and jet production phase, 
including all the potentially relevant physics like neutrino 
transport and magnetic fields. Only then, for example, will 
the question be answered under which circumstances the 
neutrino-driven wind from the transiently stable hypermas- 
sive merger remnant represents an unsurmountable obstacle 
for GRB jets. 

• The complex sequence of physical processes from the 
energy release by the merger remnant to the stage of GRB 
production involves large uncertainties, expressed by the 
product of efficiency factors /i to fi and jet collimation 
factor fn, which appears in Eq. Q. Each of these factors 
can contain still unknown dependences on the parameters of 
the NS+NS /BH binary and of its relic BH-torus system. Im- 
proved simulations of jet formation by energy release around 
post-merger BH-torus systems are necessary to clarify these 
dependences. 

• Since the compact proto-BHs found in our NS+NS 
merger simulations have significant angular momentum 
(a rcm is between 0.64 and 0.91 for the considered mergers 
of irrotational binaries; Tabled, the energy release during 
the accretion is likely to be boosted by Kerr effects. Not 
only neutrino emission and neutrino-antineutrino annihila- 
tion may play a role, but also energy release through mag- 
netic processes, e.g., by the Bland ford-Z najek (1977) mech- 
anism (see, for example. iLee et al]l2005l) . 

Much more theoretical work and numerical modeling are 
therefore needed before truly meaningful and reliable the- 
oretical predictions of GRB properties and their link to 
NS+NS/BH merger parameters will become possible. Such 
predictions need to invoke, in particular, detailed models 
of the structure of the collimated ultrarelativistic outflows 
from the merger remnants and assumptions about the pro- 
duction of the GRB emission at large distances. Detailed 
consideration of the relevant physics is clearly beyond the 
scope of this paper. A first, still very crude attempt to em- 
ploy such calculations of jets from BH-torus systems for pre- 
dicting the distributions of observable GRB p roperties has 
been attempted recently bv lJanka et alJ (|2005). 

6 SUMMARY AND CONCLUSIONS 

We have shown that nonradial oscillations and triaxial defor- 
mation and the associated tidal torques can mediate efficient 
angular momentum transfer in the compact remnant within 
the first milliseconds after the merging of two NSs. Applying 
different criteria, we estimated the mass of the torus which 
will survive the future collapse of the central, dense object 
to a BH. For the considered irrotation al systems and the e m- 
ployed relatively stiff nuclear EoS of IShen et alJ <ll998al lbl). 
we determined lower disc mass limits of 1 - 2% of the (bary- 
onic) system mass for a NS-NS (gravitational) mass ratio 
around q = 1, increasing to about 9% for systems with q 



around 0.55. Though the disc mass is more sensitive to the 
parameter q, we also observed a mild inverse relation be- 
tween torus mass and binary system mass. For the same 
value of q, slightly higher torus masses were found when 
Mjum = Mi + A/2 was lower (Fig.[2J. Non-zero temperature 
effects turned out to be important and to reduce the disc 
mass by several 10% compared to the values for the T = 
case. 

The disc masses obtained in our simulations confirm 
the viability of post-merging BH-torus systems as central 
engines of short GRBs, and they support the theoretical pos- 
sibility that thermal energy deposition above the poles of the 
BH by vv annihilation could be the primary energy source of 
short GRBs. The ultrarelativistic outflow of matter powered 
by this energy deposition was found in hydrodynamic sim- 
ulations to be highly collimated in jets with semi-o pening 
angles between about 10° and 15° iIaIov et alj|2005ft . This 
prediction seems to be supported by the observations of two 
of the four recen tly detected well-localized, short hard GRBs 
jFox et alJl2005h . 

Using typical parameters that describe the multi-step 
process from the energy release of the central engine to 
the observed GRB in this scenario, we derived estimates 
of the accreted masses and mass accretion rates (see Ta- 
ble^ from the measured GRB energies and durations, us- 
ing Eqs. Q and J5J|. For all cases we found values for M acc 
and M acc that are in the ballpark of the expectations from 
our NS+NS merger models and from models of BH accre- 
tion in post-merging BH-torus systems. Employing the same 
assumptions about the GRB engine, the extraordinarily low- 
energetic and short-duration GRB 050509b can be used to 
set strict limits for the lifetime of the post-merger hyper- 
massive neutron star. It must have collapsed to a BH in 
much less than 100ms with only little matter (~ 0.01 M@) 
remaining in a torus around the BH to partially produce 
GRB- viable outflow. This favors a symmetric or nearly sym- 
metric NS+NS binary as progenitor system, with a mass 
typical of the known galactic double neutron stars. 

Knowledge of the binary system parameters, which 
could be obtained from measurements of the gravitational 
wave chirp signal during the inspiral, might allow one to 
deduce interesting information about the properties of the 
neutron star equation of state from future short GRBs with 
determined distances. Constraining the EoS of NS matter 
by lifetime arguments for the supramassive or hypermas- 
sive remnant of NS+NS mergers, however, will require de- 
tailed relativistic models of the NS merging and of the post- 
merging evolution for different non-ze ro temperature EoSs , 
including the effects of magnetic fields flShibata et alfe oOS*). 
GRB 050509b provides particularly strict limits in such an 
analysis, but may be a rare event where an exceptionally 
low-energy GRB with very short duration could be well lo- 
calized. 
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